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We present a measurement of the Cabibbo-Kobayashi-Maskawa matrix element \V u b\, based on 
253 fb _1 of data collected by the Belle detector at the KEKB e + e~ asymmetric collider. Events 
are tagged by fully reconstructing one of the B mesons, produced in pairs from T(4S). The signal 
for b — » u semileptonic decay is distinguished from the b — > c background using the hadronic 
mass Mx, the leptonic invariant mass squared q 2 and the variable P+ = Ex — \px\- The results 
are obtained for events with the prompt-lepton momentum, p\ > 1 GeV/c, in three kinematic 
regions (1) M x < 1.7 GeV/c 2 , (2) M x < 1.7 GeV/c 2 combined with q 2 > 8 GeV 2 /c 2 , and by 
(3) P+ < 0.66 GeV/c, allowing for a comparison of the three methods. The matrix element \V u b\ 
is found to be (4.09 ± 0.19 ± 0.20 1$' is + 0.18) x 10~ 3 , where the errors are statistical, systematic 
including Monte Carlo modeling, theoretical and from shape function parameter determination, 
respectively. 

PACS numbers: 12.15.Hh,11.30.Er,13.25.Hw 



Aii accurate knowledge of the Cabibbo-Kobayashi- 
Maskawa matrix element | V u b | is crucial to test Standard 
Model predictions for CP violation. Currently, the best 
precision may be achieved by measuring the inclusive rate 
Ar M ^„(A$) of B — > X u lv decays in a restricted region of 
the phase space (A<£>) where the dominant charm back- 
ground is suppressed and theoretical uncertainties are 
minimized. The theoretical factor i?(A$) directly relates 
the inclusive rate to \V u b\, with no extrapolation to the 
full phase space: |K&| 2 = Ar uf „(A$)/i?(A$). Here we 
report measurements of Ar tl £ l/ (A$) for several choices of 
A$ and derive the corresponding values of \V u b\. 

The measurements are made with a sample of events 
where the hadronic decay mode of the tagging side B me- 
son, Stag, is fully reconstructed, while the semileptonic 
decay of the signal side B meson, i? s i g , is identified by the 
presence of a high momentum electron or muon. B de- 
notes both charged and neutral B mesons. This method 
allows the construction of the invariant masses of the 
hadronic (Mx) and leptonic (y/<^) system in the semilep- 
tonic decay, and the variable J+ = Ex — \px \ , where 
Ex is the energy and \px\ the magnitude of the three- 
momentum of the hadronic system. These inclusive kine- 
matic variables can be used to separate the B — > X u tv 
decays from the much more abundant B — > XJtv decays. 
Three competing kinematic regions (A$) were proposed 



by theoretical studies [1, 2], based on the three kinematic 
variables, and are directly compared by this analysis. 
The value of \V u b\ is extracted using recent theoretical 
calculations [2, 3] that include all the currently known 
contributions. Mx and q 2 selections were already used 
to to extract \V u b\ [4, 5]. The present analysis is the first 
one to use P+ and to directly compare the three methods. 

The data were collected with the Belle detector [6] at 
the asymmetric-energy KEKB storage ring [7]. The re- 
sults presented in this paper are based on a 253 fb _1 
sample recorded at the T(4S) resonance, which con- 
tains 275 x 10 6 B~B pairs. An additional 28 fb" 1 sam- 
ple taken at a center-of-mass energy 60MeV below the 
T(4S f ) resonance is used to subtract the background from 
e + e~ — > qq (q = u, d, s, c). 

Monte Carlo (MC) simulated events were used to de- 
termine efficiencies as well as signal and background 
distributions. The detector simulation was based on 
GEANT [8]. To model B -> X u iv we use the Evt- 
Gen generator [9] with various models, where X u is -k 
or p [10], an excited X u state [11], or a non-resonant 
multiparticle final state [12]. The B — > X c lv transitions 
are simulated according to the QQ generator [13]. For 
the two dominant contributions, D*£v and Dtv, we use 
an HQET-based parametrization of form factors [14] and 
ISGW2 model [11], respectively. For the D** we use 
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ISGW2 model and for sub-components D\ and D\ set 

B ^ D bZd^v IIV = °- 35 ± °- 23 - Thc motion of the b 
quark inside the B meson is implemented with the intro- 
duction of a shape function [12, 15] that describes the b 
quark momentum distribution inside the B meson. 

The i?tag candidates are reconstructed in the modes 
B -> D^n/p/ai/D^, D° -> K+ir', K+tt-j: , 
K+tt+it-tt-, KgTT , iT^TT+Tr", Kg7T^7T~7T® and K + K~ , 
D~ -> K+tt-tt-, K+n-ir-ir , K^tt~, K^-ir , 
KgTr~ir~Tr + and A' + A'~7r~, and Df — ► K° S K + and 
if + A' _ 7r + . D* mesons are reconstructed by combining 
a D candidate and a soft pion or photon. (Inclusion of 
charge conjugate decays is implied throughout this pa- 
per.) The selection of B tag candidates is based on the 
beam-constrained mass, M bc = ^/E^^/c 4 — p* B 2 /c 2 , 
and the energy difference, AE = Eg — ££ eam . Here 
^beam = \T S I^ — 5.290 GeV is the beam energy in the 
e + e~ centcr-of-mass system (cms), and p* B and E* B are 
the cms momentum and energy of the reconstructed B 
meson. (Throughout this paper the variables calculated 
in the cms are denoted with an asterisk.) 

The combinatorial background from jet-like e + e~ — > 
qq processes is suppressed by an event topology re- 
quirement based on the normalized second Fox- Wolfram 
moment i?2 < 0.5 [16], and for some modes also by 
| cos6>* hrust | < 0.8, where #* hrust i s the angle between the 
thrust axis of thc Stag candidate and that of the rest 
of the event. To minimize the fraction of events with 
incorrect separation of tag and signal sides while main- 
taining high signal efficiency, a loose selection require- 
ment of M bc > 5.22 GeV/c 2 and -0.2 < AE < 0.05 GeV 
is made. If an event has multiple Stag candidates, we 
choose the one having the smallest x 2 based on AE, the 
D candidate mass, and the D* — D mass difference if 
applicable. 

For events tagged by fully reconstructed i? ta g candi- 
dates, we search for electrons or muons from semileptonic 
decays of B slg . We require a lepton with momentum p* t 
exceeding 1 GeV/c in thc laboratory polar angular region 
of 26° < 9 < 140°. Lcptons from J/ip decay, photon con- 
version in the material of the detector, and ir° decay are 
rejected based on the invariant mass they form in com- 
bination with an oppositely charged lepton and for elec- 
tron candidates also with an additional photon. When 
the Btag candidate is charged, we also require the lepton 
charge to be consistent with that from prompt semilep- 
tonic decay The signal yield is obtained by fitting the 
Mb c distribution to the sum of an empirical parametriza- 
tion of the combinatorial background shape [17] plus a 
signal shape [18] that peaks at the B mass and taking 
the part of the signal that lies in thc "signal region," 
M bc > 5.27GcV/c 2 , as shown in Fig. 1(a). Thc cut- 
off for Mbc reduces the uncertainty from the incorrect 
assignment of tag and signal sides in signal events. 

The B — > X u lv signal events are selected by re- 
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FIG. 1: (a) Distribution in A/bc (data) of Stag candidates 
in events satisfying _B s i g selection, (b) Mx distribution for 
events with q 2 > 8 GeV 2 /c 2 , with fitted contributions of 
B -> Xclv and B -» X u tv. 

moving poorly measured soft charged tracks and impos- 
ing several additional requirements to reject poorly re- 
constructed events and suppress the B — > X c lv back- 
ground. We require that the event contain exactly one 
lepton and have zero net charge and that the invariant 
mass squared of the missing four-momentum m^-^ = 

(pt(4S) ~ Ps t a g ~ Px - Pt) 2 (PT(4S), Ps tag and Px are 
four-momenta of the T(4S'), B tag , and hadronic system 
PO, respectively) be within -1 < m 2 liss < 0.5GeV 2 /c 4 . 
To suppress the B — > X c £v background, events with 
a K ± or Kg candidate on thc signal side are rejected 
(kaon veto). To reject events containing a X£, we re- 
quire that the angle between the missing momentum 
and the direction of any candidate, reconstructed 
in the K\ detector, be greater than 37 degrees. We 
also reject B° — » D* + (~D events by detecting the slow 
pion (ir s ) from D* + — > D°tt+ and deducing from its 
momentum the momentum of the D* + . The missing 
mass squared , D *-> = (pb — Pd* — Pi) 2 is calcu- 

lated from the reconstructed quantities, and events with 
to 2 . , n „N > — 3GeV 2 /c 4 are rejected. 

miss \D * ) I J 

Finally, the kinematic variables Mx and P+ are calcu- 
lated from the measured momenta of all charged tracks 
and energy deposits of all neutral clusters in the elec- 
tromagnetic calorimeter that are not used in the Stag 
reconstruction or for the lepton candidate. The four- 
momentum of thc leptonic system is calculated as q = 
Pt(4S) ~ P-Btag ~ Px- The distributions of events in Mx 
and P+ are obtained by fitting the Mb c distribution, as 
described above, in bins of Mx and P+. Figures 1(b), 
2(a) and 3(a) show the resulting Mx and P+ distribu- 
tions. We define three kinematic signal regions (A$) for 
events where the prompt lepton has p\ > 1 GeV/c: P + < 
0.66 GeV/c, M x < 1.7 GeV/c 2 , and M x < 1.7 GeV/c 2 
combined with q 2 > 8 GeV 2 /c 2 . These three regions 
are denoted as P+, Mx and Mx/q 2 , respectively. To 
minimize the systematic effects of uncertainties in lepton 
selection and full reconstruction, we normalize the par- 
tial rate for each signal region with thc total semileptonic 
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TABLE I: JVJ! 
signal regions. 



-scT U ! F an d r b-^u f° r the three kinematic 
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To extract the raw number of signal events, N[™ u , we 
fit the Mx and P+ distributions with MC-determined 
shapes for B — > X u £v and P — > ATc^ and subtract the 
P — ► X c ^z/ contribution. The results for the Mx/q 2 and 
P+ regions are shown in Figs. 1(b) and 2(a), respectively. 
MC simulation is used to estimate the conversion factor 
F of the observed number of events N£™ u to the num- 
ber of signal events produced in the region in question 
and observed anywhere, and to estimate the efficiency 
for these events, e^ u . 

N a \ = (9.14±0.05) x 10 4 is the number of events having 
at least one lepton with p\ > 1 GeV/c, determined from a 
fit to the corresponding Mb c distribution (Fig. 1(a)), and 
corrected for the expected fraction of background events 
from non-semileptonic decays (14.0%), as estimated by 
MC simulation. The factor ef' ec /ef r "^ c " accounts for a 
possible difference in the P tag reconstruction efficiency 
in the presence of a semileptonic or B — > X u £v decay; 
£ f/ £ e~~* u i s t ne ratio of both lepton identification effi- 
ciencies and fractions of semileptonic decay leptons with 
P*e > lGcV/c, in the whole kinematic phase space for 
semileptonic decays, and within the kinematic signal re- 
gion for signal events. The product of efficiency ratios 



ulation. 



"frcc/°frcc 



X £' 



l /4 



is obtained from MC sim- 



Table I summarizes the results for 



„fa — >i 

-se\ 



F and rl_^ u for all three signal regions, where the error 
in N£™ u is statistical only. Inserting these values in Eq. 
1, we obtain the three values of W. As both numerator 
and denominator of W have been obtained from the same 
tag sample, the B°/B + weightings are the same, and W 
has no dependence on lifetimes. Multiplying W by the 
average measured semileptonic rate T(X£v) = B(B — > 
X£v)/t b , obtained from B(B -> Xiv) = 0.1073 ± 0.0028 
and tb = (1.604 ± 0.016) ps [19], gives the average 
Ar u £„(A$). The results with relative errors are given 
in Table II. 

We divide the experimental error into four categories: 
statistical, systematic, b — > c and b — > u MC mod- 
eling errors, and summarize them in Table II for the 
three Ar„^ l/ (A$) measurements. The two modeling er- 
rors include the uncertainty in signal event extraction, 
efficiency and unfolding factor determination due to the 
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FIG. 2: (a) The P+ distribution for the selected events, with 
fitted contributions from B — > X c £v and B — > X u tv, (b) 
P+ distribution (symbols with error bars) after subtracting 
B — » X c lv, with fitted B — > X u £v contribution (histogram). 



TABLE II: Partial rates to the three kinematic signal regions 
with relative errors (in %). 





xr utv {A$) 


stat 


syst 


b — > u 


b -» c 


Mx/q' 


5.24 x 10" 4 ps" 1 


10.0 


8.9 


6.2 


5.3 


M x 


7.71 x 10~ 4 ps" 1 


9.1 


7.1 


6.1 


2.2 


P+ 


6.89 x 10 -4 ps" 1 


9.4 


9.3 


6.4 


8.7 



choice of specific theoretical models and values of the 
parameters used in our MC predictions. For signal 
B — > X u lv MC, the shape function parameters A SF = 
(0.66 ± 0.15) GcV/c 2 and Xf F = -(0.40 ±0.20) GcV 2 /c 2 
were varied within the stated limits, taking into account 
the negative correlation between them [20] . To take into 
account the uncertainty of the prediction in Ref. [12], we 
use a factor of two larger error for A SF than was deter- 
mined in Ref. [20]. For B — > X c £v MC, the uncertainty 
due to our limited knowledge of branching fractions is 
studied by varying the contributions of Dlv and D*£v 
and the relative fraction of narrow states D\ and Pj that 
contribute to D**£v to estimate the modeling error of the 




M s (GeV/c") 



M x (GeV/c") 



FIG. 3: Mx distribution (no q 2 requirement) with fitted con- 
tributions from X c £v and X u £v: (a) before and (b) after 
subtracting the X c £v contribution (symbols with error bars), 
shown with the prediction for X u £v (MC, histogram). 
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D** region. The uncertainty from form factor modeling 
in Diu and D*£v was studied by varying the parameters 
p 2 D = 1.15 ± 0.16 and p\ = 1.56 ± 0.13 within their er- 
rors [19]. The validity of the B — > X c lv simulation was 
tested on a B — ► X c £i/ enhanced control sample, where 
all selection requirements are applied but with the kaon 
veto reversed. The kinematic distributions of this con- 
trol sample are accurately described by the simulation. 
Other sources of uncertainties, namely limited MC statis- 
tics, extraction of rf_^ u , fitting procedure and imperfect 
detector simulation are combined in the systematic error. 
The uncertainties due to inaccurate simulation of track- 
ing, particle identification, and cluster finding are esti- 
mated by varying for each source the efficiency within 
the expected error and taking the maximum change in 
Ar u ^„(A$) as the error. For each of these sources the 
effects on simulated b — > u and b — > c events are cor- 
related, and the associated shifts are summed linearly. 
The net contributions from the three sources are then 
summed in quadrature. 

The CKM matrix element \V u b\ is obtained directly 
from the partial rate using \V ub \ 2 = Ar„^(A$)/i?(A$). 
i?(A$) is the theoretical prediction of AT U £^(A^>), the 
partial rate with a prompt lepton with p* e > 1 GeV/c, di- 
vided by |Vub| 2 - The values of R (in ps _1 ) are calculated 
to be 23.7 ± 2.0(SP)±|;|(th.), 46.1 ± 4.2(SF)Tj'j(th.) 
and 39.4 ± 4.5(SF)^-8(th.) for the M x /q 2 , M x and 
P + signal regions, respectively The i?(A$) values and 
their errors (SF) arc calculated using the shape function 
scheme [15] parameters m b (SF) = (4.60 ±0.04) GcV/c 2 
and a4(SF) = (0.20 ±0.04) GeV 2 /c 2 with correlation 
coefficient p = —0.26, obtained from the result of a 
global fit to moments of both b — > civ and b — ► S7 
distributions [21]. While the dependence of i?(A$) on 
/i 2 (SF) is small, we can approximate the dependence on 
m & (SF) as R/R{m° b ) - 1 = fc(A$) • (m b /m^ - 1), where 
ml = 4.60 GcV/c 2 and fc(A$) is found to be 2.09, 2.29 
and 3.00 for the Mx/q 2 , Mx and P + signal regions, re- 
spectively. The theoretical error of R (th.) is estimated 
by varying the sublcading shape functions (four models), 
the matching scales fi^, (Uj, /2 and weak annihilation [15]. 
The values of \V u b\ with errors are given in Table III. The 
total error on \V ub \ is 10%, 9% and 11% for Mx/q 2 , M x 
and P+ regions, respectively. When the shape function 
parameters and R are better determined, \V ub \ can be 
recalculated from Ar tl ^ l/ (A$) shown in Table II. 

The precision of the \ V ub \ determination is better than 
previous measurements [4, 5, 22], owing to the use of 
larger data sample, better shape function parameter de- 
termination and improved theoretical predictions [2, 3]. 
We find that the usage of the variable P + is more sensi- 
tive to b — ► c modeling and shape function parametriza- 
tion than the other two methods and will become com- 
petitive in the future when the theoretical error of R 
dominates. No significant experimental nor theoretical 
improvement was observed by applying the additional 



TABLE III: Values for \V ub \ with relative errors (in %) for the 
three kinematic signal regions. Shape function parameters 
used in the calculation are m b (SF) = (4.60 ± 0.04) GeV/c 2 
and a4(SF) = (0.20 ± 0.04) GeV 2 /c 2 . 





\V ub \x 10 3 


stat 


syst 


b — > u 


b -> c 


SF 


th. 


Mx/q' 


4.70 


5.0 


4.4 


3.1 


2.7 


4.2 


+4.8 

-5.2 


M x 


4.09 


4.6 


3.5 


3.1 


1.1 


4.5 


+3.5 
-3.8 


P+ 


4.19 


4.7 


4.6 


3.2 


4.4 


5.8 


+3.4 

-3.5 



selection q 2 > 8 GeV 2 /c 2 to the Mx analysis. Tak- 
ing correlations into account, we find that the difference 
between \V ub \ values for Mx/q 2 and Mx regions has a 
significance of 2.7er. We conclude that the results are 
consistent within errors, but we do not rule out possible 
effects of duality violation or weak annihilation contri- 
bution. We chose the Mx signal region result for our 
|V^b| determination, since it includes the largest portion 
of phase space and is least affected by the uncertainties: 
\V ub \ =(4.09 ± 0.19 ± 0.20 to.lt ± 0.18) x 10" 3 , where 
the errors are statistical, systematic with MC modeling, 
theoretical and from shape function parameter determi- 
nation, respectively. The effectiveness of \V ub \ measure- 
ments using full reconstruction tagging is clear (Figs. 2(b) 
and 3(b)). 
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